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The question of Fermi liquid vs. non-Fermi liquid behavior induced by 
strong correlations is one of the prominent problems in metallic local mo- 
ment systems. As standard models for such systems, the SU(N)xSU(M) 
Anderson impurity models exhibit both Fermi liquid and non-Fermi liquid 
behavior, depending on their symmetry. Taking the Anderson model as 
an example, these lectures first give an introduction to the auxiliary boson 
method to describe correlated systems governed by a strong, short-range 
electronic repulsion. It is then shown how to include the relevant low-lying 
excitations (coherent spin flip and charge fluctuation processes) , while pre- 
serving the local gauge symmetry of the model. This amounts to a conserv- 
ing T-matrix approximation (CTMA). We prove a cancellation theorem 
showing that the CTMA incorporates all leading and subleading infrared 
singularities at any given order in a self-consistent loop expansion of the 
free energy. As a result, the CTMA recovers the correct infrared behavior 
of the auxiliary particle propagators, indicating that it correctly describes 
both the Fermi and the non Fermi regimes of the Anderson model. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.28.+d, 75.20.Hr 

1. Introduction 

It is a remarkable feature of interacting, itinerant fermion systems that 
at low temperatures T they behave in general in much the same way as a 
noninteracting Fermi gas, even though the interaction may be strong. An 
extremely successful description of this phenomenon, known as Fermi liq- 
uid (FL) behavior, is provided by the notion of quasiparticles, which was 
established by Landau's phenomenological Fermi liquid theory [1]. The 
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key assumption is that, as the interaction is continuously turned on, there 
exists a 1:1 correspondence between the low energy eigenstates of the inter- 
acting system and the single-particle states of the free Fermi gas. Therefore, 
the low-lying interacting states may be described approximately as single- 
particle states or quasiparticles, whose decay rate 1/r is small compared to 
their excitation energy u, 1/r <C u, and which are characterized by the same 
quantum numbers as the noninteracting states. As a consequence, Fermi 
liquids exhibit the same low-T thermodynamics as a noninteracting Fermi 
system, e.g. a linear in T specific heat c = 7T and a constant Pauli param- 
agnetic susceptibility Xo- However, the effective mass and other parameters 
may be renormalized by the interaction, resulting in an enhancement of the 
specific heat coefficient 7 and the susceptibility \o- It is & t the heart of the 
quasiparticle picture that at low T the Pauli exclusion principle substan- 
tially reduces the phase space available to quasiparticle scattering. This 
blocking mechanism is effective as long as the quasiparticle interaction is 
shortranged in space and time, which is usually the case in three dimen- 
sions because of screening. It also implies that the quasiparticle scattering 
rate vanishes as 1/r oc (uj 2 + T 2 ) in the limit u>, T — > 0, thus providing a 
microscopic justification for the basic assumption of FL theory and leading 
to an interaction contribution to the electrical resistivity which behaves as 
Ap oc T 2 . Obviously, the Pauli principle as the origin of FL behavior is very 
robust, which explains the almost ubiquitous presence of a FL ground state 
in interacting Fermi systems and the broad success of Fermi liquid theory. 

In this light it is all the more exciting that in recent years a number of 
new alloys have been discovered which exhibit striking deviations from this 
usual behavior. These systems have in common that a localized, degenerate 
degree of freedom, the magnetic moment of a magnetic ion, is dynamically 
coupled to a continuum of conduction electron states. In general, such a 
coupling generates the Kondo effect, characterized by resonant spin flip scat- 
tering of electrons at the Fermi surface off the local moment. Concomitantly, 
the conduction electron spin flip rate initially increases logarithmically as 
the temperature is lowered, passes through a maximum at a characteristic 
scale, the Kondo temperature Tk, and approaches zero as T — > 0, because 
the effective local moment becomes screened by the conduction electron 
spins. Thus, even for many strongly correlated systems of this type a Fermi 
liquid description applies below Tx, with usually a strongly enhanced quasi- 
particle effective mass, lending the term "heavy fermion systems" to these 
materials. 

Completely new physics may arise, however, if the quenching of the local 
moments is inhibited. Two different mechanisms for the appearance of a 
non-Fermi liquid ground state have been put forward: 

(1) Proximity of a quantum phase transition (QPT) to an antiferromag- 
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netically ordered state [2]-[5] as a function of a dopant concentration 
x or of pressure. Near the QPT the quantum critical fluctuations be- 
come longranged in space and time and can, thus, mediate a longrange 
quasiparticle interaction, leading to a breakdown of FL theory. There 
are indications for this spatially extended mechanism to be realized 
near the QPT of certain Ce based compounds like CeCu§- x hxi x [6, 7], 
CeCu 2 Si 2 [8], and CePd 2 Si 2 [9, 10]. 

(2) Two-channel Kondo effect (2CK) [11, 12]. The local magnetic mo- 
ment is coupled to two exactly degenerate conduction electron chan- 
nels. Because of a frustration effect between the screening of the local 
moment by the different conduction channels the moment quenching 
cannot be complete, leading to a nonvanishing conduction electron 
spin scattering rate even at the lowest temperatures and subsequently 
to a breakdown of FL behavior. It has been suggested [13] that this 
mechanism, based on single-ion dynamics rather than longrange fluc- 
tuations, may be realized predominantly in U based materials with 
cubic symmetry about the magnetic ion, such as Yi_ x U x Pd3 [14] or 
UCu5_ x Pt x [15], which do not exhibit a QPT. 

In both scenarios the wealth of experimental data showing non-FL be- 
havior at low temperatures is not consistently explained by the present 
theories. Open questions within the QPT picture include, e.g., whether the 
local impurity dynamics competing with the magnetic ordering can play 
a role, and how the transition from the spin screened heavy FL phase to 
the magnetically ordered phase occurs alltogether. In the 2CK mechanism, 
on the other hand, inter-impurity interactions could modify the single ion 
behavior. Exact solution methods as well as numerical simulations have 
provided important progress in our understanding of strongly correlated 
quantum impurity systems. However, their applicability is essentially re- 
stricted to problems involving only a single impurity, owing to integrability 
conditions or limitations in the numerical effort, respectively. Therefore, 
more generally applicable theoretical techniques are called for. 

In the present work we focus on the single-ion dynamics. We develop 
a standard field theoretical method, based on an auxiliary particle or slave 
boson representation, which describes the quantum impurity dynamics in a 
controlled way and at the same time has the potential of being extended to 
problems of many impurities on a lattice. As a standard model of strongly 
correlated electrons which, depending on its symmetry, exhibits both FL and 
non-FL behavior, we consider the SU(N) xSU(M) Anderson impurity model 
of a local, iV-fold degenerate degree of freedom, coupled to M identical 
conduction bands. 
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In order to set the stage for the more formal development of the the- 
ory, in the following section we will briefly review the striking differences 
in the phenomenology of the single-channel and the multi-channel Kondo 
effects. In section 3 the slave boson representation is introduced, which 
provides a particularly compact formulation of the SU(N)xSU(M) Ander- 
son model. We also discuss why the presence of FL or non-FL behavior 
in a given quantum impurity system can already be seen from the singular 
infrared dynamics of the auxiliary particles. Section 4 contains a critical as- 
sessment of earlier approximate slave boson treatments. This will motivate 
our conserving slave boson approach (conserving T-matrix approximation, 
CTMA), which is developed in section 5. As will be seen, the results pro- 
duced from this theory are in very good agreement with known exact proper- 
ties of the model. Conclusions are drawn in section 6. In the appendices we 
prove a cancellation theorem for non-CTMA diagrams, which justifies the 
CTMA on formal grounds, and derive in detail the self-consistent CTMA 
equations. 

2. Single and multi channel Kondo effect and possible physical 

realizations 

In this section we briefly discuss how the single-channel and the two- 
channel Kondo effects may arise in magnetic metals, if the interaction be- 
tween the local moments can be neglected. We first discuss the usual mag- 
netic, single-channel Kondo effect. A local moment is generated by an 
atomic f or d level whose energy Ed lies far below the Fermi energy ep = 
and whose electron occupation number is effectively restricted to rid < 1 by 
a strong Coulomb repulsion U between two electrons in the same orbital. 
While the angular momentum degeneracy of the level is usually lifted by 
crystal field splitting, in the absence of a magnetic field a twofold degeneracy 
of a level occupied by one electron is guaranteed by time reversal symmetry 
(Kramers doublet), corresponding to the spin quantum numbers m = ±1/2 
of the electron. In addition, there is a hybridization matrix element V be- 
tween the atomic orbital and the conduction electron states. Such a system 
is described by the single impurity Anderson Hamiltonian 

H = E £ k c l c k* + E * E d Ua + v E(4> + /l - c - ) + ud \ d A d v ( 2J ) 

k,a a k,a 

where ct and dl are the creation operators of a conduction electron with 

ka 

dispersion and of an electron in the local orbital with spin a, respectively. 
The low energy physics of this system is dominated by processes of second or- 
der in V, by which an electron hybridizes with the conduction band and the 
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impurity level is subsequently filled by another electron, thereby effectively 
flipping the impurity spin. Thus, in the region of low excitation energies, 
the Anderson Hamiltonian (2.1) may be mapped onto the s-d exchange (or 
Kondo) model [16], the effective coupling between the impurity spin and the 
conduction electron spin always being antiferromagnetic: J=|y| 2 /|E^|>0 
(U 3> l-E^I). These models have been studied extensively by means of 
Wilson's renormalization group [17], by the Bethe ansatz method [18, 19] 
and by means of a phenomenological Fermi liquid theory [20]. In this way 
the following physical picture has emerged: The model contains a dynam- 
ically generated low temperature scale, the Kondo temperature, which is 
expressed in terms of the parameters of the Anderson Hamiltonian (2.1) 
as T K = D(NT/D)( M Wexp{-7rE d /(NT)}, with Af(0) and D = 1/Af(0) 
the density of states at the Fermi energy and the high energy band cutoff, 
respectively. T = irV 2 N(0) denotes the effective hybridization or d-level 
broadening and N, M are the degeneracy of the local level and the number 
of conduction electron channels (see below). In the intermediate tempera- 
ture regime, T ~ Tk, resonant spin flip scattering of electrons at the Fermi 
surface off the local degenerate level leads to logarithmic contributions to the 
magnetic susceptibility, the linear specific heat coefficient and the resistivity 
(p(T) = p(0) + Ap(T)), X (T), 7(T), Ap(T) oc -\n(T/T K ), and to a break- 
down of perturbation theory at T ~ Tk- Below Tk a collective many-body 
spin singlet state develops in which the impurity spin is screened by the 
conduction electron spins as lower and lower energy scales are successively 
approached, leaving the system with a pure potential scattering center. The 
spin singlet formation is sketched in Fig. 1 a) and corresponds to a vanishing 
entropy at T = 0, S(0) = 0. It also leads to saturated behavior of physical 
quantities below Tk, like xCO = const, c(T)/T = const, and Ap(T) oc T 2 , 
i.e. to Fermi liquid behavior. 

As an example of possible two-channel Kondo systems we discuss the 
uranium based compounds mentioned in the introduction. The U 4+ ions 
have nominally a 51 2 configuration, i.e. an even number of electrons, which 
does not allow for a Kramers degenerate ground state because of integer 
total spin. However, in the cubic crystal symmetry of these materials the 
orbital degeneracy may be not completely lifted, so that there can be an 
approximate twofold degeneracy of the U 4+ ground state, corresponding 
to two different orientations of the electrical quadrupole moment of the 5f 
orbital in the lattice (quadrupolar Kondo effect) [12, 13]. This degree of 
freedom may be flipped by scattering of conduction electrons (which in the 
cubic symmetry also have a twofold angular momentum degeneracy). The 
conduction electron spin is conserved in this scattering process, leaving it as 
a Kramers degenerate scattering channel degree of freedom, which we will 
label by fi = 1, M, M = 2. Describing the orbital degree of freedom as 
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Fig. 1. Sketch of the renormalization group for a) the single-channel Kondo model 
(local moment compensation) and b) the two-channel Kondo model (local moment 
over-compensation). Small arrows denote conduction electron spins 1/2, a heavy 
arrow a localized spin 1/2. The curved arrows indicate successive renormalization 
steps. 

a pseudospin 1/2, labelled by the quantum number a = 1, . . . , N, N = 2, 
in analogy to the magnetic Kondo effect, we arrive at the SU(2)xSU(M) 
symmetric Kondo model, 

H=Ye t cl cr +J Y ct S-T aa <c r (2.2) 

where 5 is the local pseudospin operator and r CTCT / the vector of Pauli ma- 
trices. To keep the naming uniform, we will refer to the orbital degree of 
freedom as the (pseudo)spin or local moment, a, in analogy to the magnetic 
Kondo effect, and to the physical electron spin as the channel degree of free- 
dom, jx. In the multi-channel case, too, the conduction electrons of each 
channel separately screen the impurity moment by multiple spin scattering 
at temperatures below the Kondo scale Tk- However, in this case, the lo- 
cal moment is over-compensated, since the impurity spin can never form 
a singlet state with both conduction electron channels at the same time in 
this way, as can be seen in Fig. lb). As a consequence of this frustra- 
tion, there is not a unique ground state, leading to a finite residual entropy 
[21, 22] at T = of 5(0) = fc^lriV^ in the two-channel model. In particular, 
the precondition of FL theory of a 1:1 correspondence between interacting 
and non-interacting states is violated. As a consequence, characteristic sin- 
gular temperature dependence [21, 23] of physical quantities persists for 
T ~ Tk down to T = 0: X (T) oc -\n(T/T K ), c(T)/T oc -ln(T/T K ) and 
p(T) — p(0) oc —\JT/Tk- Note, however, that this behavior may be changed 
by any crystal field splitting of the quadrupolar non-Kramers ground state 
doublet. 

In order to apply standard field theoretical methods to the multi-channel 
Kondo model it is convenient to consider it as the low-energy limiting case of 
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a corresponding Anderson model as discussed for the single-channel case. 
Here, in addition, the conservation of the channel degree of freedom has 
to be guaranteed. This can be implemented in an elegant way using an 
auxiliary boson representation to be discussed in the next section. 

3. Auxiliary particle representation 

As discussed above, the local level of a quantum impurity in the limit of 
infinitely strong local Coulomb repulsion U between electrons in the same 
level allows only for at most single electron occupation of the level, < 1. 
One should note that for a realistic finite value of U the low-energy physics 
of the model is effectively still confined to the part of the Hilbert space 
without multiple occupancy. Therefore, the model Eq. (2.1) in the limit 
U — ► oo is the generic model for the physics of quantum impurities at large 
U in general. 

A powerful technique for implementing the projection in Hilbert space 
caused by a large Coulomb repulsion U is the method of auxiliary particles 
(slave bosons, pseudofermions) [24]. Each Fock state \a) of the impurity 
is assigned a creation operator, which can be envisaged as creating the 
state out of a vacuum state |vac) without any impurity level at all, \a) = 
at, | vac). (E.g., for a single orbital there are four such states, |0) (empty 
orbital), | |) or | j) (orbital occupied by a single electron with spin f or 
I) and 1 2) (level occupied by two electrons with spin f and J,).) Due to 
the requirements of Fermi statistics, the creation operators a) a are Fermi 
(Bose) operators for the states holding an odd (even) number of electrons 
(or vice versa). The physical state corresponds to the sector of Fock space 
with exactly one auxiliary particle, ^2 a n a = 1, where the 
occupation number operator of particles a. Compared to alternative ways of 
effecting the projection, the auxiliary particle method has the advantage of 
making available the powerful machinery of quantum field theory, provided 
the constraint on the total auxiliary particle number can be incorporated 
in a satisfactory way. 

For the quantum impurity models of the Anderson type introduced in 
the preceding section, only particles creating empty and singly occupied 
states are needed. We define N pseudofermion creation operators /J for 
each of the singly occupied states (labelled by a = 1,2, ... ,N) and M 
boson creation operators for each of the empty states created when an 
electron hops from the impurity into the ^-th conduction electron band 
(labelled by /x = 1, 2, . . . , M). In terms of these operators the Hamiltonian 



8 



CTMA PRINTED ON FEBRUARY 1, 2008 



of the SU(N)xSU(M) Anderson model Eq. (2.1) takes the form 

H = J2 ^IAo + E * E &f' + V J2 ^\J»U + h.c.) (3.1) 

In order for H to be SU(M) invariant, the slave boson multiplet bp trans- 
forms according to the conjugate representation of the SU(M). In addition, 
the operator constraint 

Q = Y,ftu + Y, h lh = l (3-2) 

has to be satisfied at all times. One might interpret the constraint as a 
statement of charge quantization, with the integer Q the conserved, quan- 
tized charge. Similar to quantum field theories with conserved charges, the 
charge conservation is intimately related to the existence of a local gauge 
symmetry. Indeed, the system defined by the Hamiltonian Eq. (3.1) is in- 
variant under simultaneous local U(l) gauge transformations j c — > f a e l ^ T \ 
bp — ► bpe 1 ^^ , with <p(r) an arbitrary time dependent phase. 



3.1. Exact projection onto the physical Hilbert space 



While the gauge symmetry guarantees the conservation of the quantized 
charge Q, it does not single out any particular Q, such as Q = 1. In order 
to effect the projection onto the sector of Fock space with Q = 1, one may 
use a procedure first proposed by Abrikosov [25]: Consider first the grand- 
canonical ensemble with respect to Q, defined by the statistical operator 

PG = ^^ H+XQ \ (3.3) 

where Zq = tr[exp{— (3{H + AQ)}] is the grand-canonical partition function 
with respect to Q, —A is the associated chemical potential, and the trace 
extends over the complete Fock space, including summation over Q. The 
expectation value of an observable A in the grand-canonical ensemble is 
given by 

(A) G = ti[p G A]. (3.4) 

The physical expectation value of A, (A), is to be evaluated in the canonical 
ensemble where Q = 1. It can be calculated from the grand-canonical 
ensemble by differentiating with respect to the fugacity £ = e _/3A and taking 
A to infinity [26], 



£ tr[Ae-^ H+x ®} (QA)g 
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Projecting operators acting on the impurity states: — We list two important 
results, which follow straightforwardly from Eq. (3.5): First, the canonical 
partition function in the subspace Q = 1 is 

Z C = lim tr[Qe-^ H+x ^-^] 

A— >oo 

= lim (e l3X (Q) G (X))Z Q=0 , (3.6) 

A— >oo 

where the subscripts G and C denote the grand-canonical and the canonical 
(Q = 1) expectation value, respectively. Second, the canonical Q = 1 
expectation value of any operator A which has a zero expectation value in 
the Q = subspace, A\Q = 0) =0, is given by, 

{A)c = lim (3.7) 

Note that A\Q = 0) =0 holds true for any physically observable opera- 
tor acting on the impurity. Examples are the physical electron operator 
d-lur = fa^p. or the local spin operator S = J2aa> ^fl^cra' fa' ■ I n this case 
the operator Q appearing in the numerator of Eq. (3.5) is not necessary to 
project away the Q = sector. In particular, the constrained ci-electron 
Green's function is given in terms of the grand-canonical one (Gd(to,T, A)) 
as 

G d {u) = lim — (3.8) 

In the enlarged Hilbert space (Q = 0, 1, 2, ...) Gd(to, T, A) may be expressed 
in terms of the grand-canonical pseudo-fermion and slave boson Green's 
functions using Wick's theorem. These auxiliary particle Green's functions, 
which constitute the basic building blocks of the theory, are defined in imag- 
inary time representation as 

9f*(n-T 2 ) =-{T{f a {r l )fl{r 2 )}) G (3.9a) 

Gbpin - r 2 ) = -<T{6p(ri)6t(7a)}) G , (3.9b) 

where T is the time ordering operator. The Fourier transforms of Qf^ may 
be expressed in terms of the exact self-energies E as 

g f , b (iu n ) = {[Ql^n)]- 1 - S /)6 (^ n )} _1 (3.10) 

where 

g%{iu n ) = (icon -E d - A)" 1 (3.11a) 



g%{iu n ) = (ico n ~ A)" 1 (3.11b) 
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Since as a consequence of the projection procedure A — > oo the energy 
eigenvalues of H + \Q scale to infinity as XQ, it is useful to shift the zero 
of the auxiliary particle frequency scale by A (in the Q = 1 sector) and to 
define the "projected" Green's functions as 

G f>b (u) = lim Gf, b (u + X) (3.12) 

A— *oo 

Note that this does not affect the energy scale of physical quantities (like 
the local d electron Green's function), which is the difference between the 
the pseudo-fermion and the slave-boson energy. 

Canonical expectation values of conduction electron operators: — The 
canonical (i.e. projected onto the Q = 1 subspace), local conduction electron 
Green function is given as 

G CfMa (iuj n ) = {[Gc^iiuJn)]' 1 - E c/10 -(iu; n )j (3.13) 

with 

G%A^ n ) =Y.G%Aki^n) =E(^ + ^- e fe) _1 . ( 3 - 14 ) 

k k 

where fi c is the chemical potential of the conduction electrons. The canon- 
ical, local c-electron self-energy, £ c/KT (zu; n ), cannot be obtained from the 
grand-canonical one by simply taking the limit A — > oo, since the c-electron 
density has a non-vanishing expectation value in the Q = subspace. How- 
ever, it follows immediately from the Anderson Hamiltonian (2.1), that the 
exact, canonical conduction electron t-matrix t a ^(iuj), defined by G c = 
G®[l + tG®], is is proportional to the full, projected <i-electron Green's func- 
tion, t^ a (iuj) = iV^Gdfxa- Thus, we have as an exact relation, 

G c „ a (iu n ) = G° CfUT (ioj n ) [l + \V\ 2 G dlMT (iu n )G° c ^(iLu n )] , (3.15) 

and by comparison with Eq. (3.13) we obtain the local conduction electron 
self-energy respecting the constrained dynamics in the impurity orbital, 

T.c^ilUn) = | T/ | 2r , f . 77; T. 7 • (3.16) 

Using phenomenological Fermi liquid theory [20] and also by means of per- 
turbation theory to infinite order in the on-site repulsion U [27] it has 
been shown for the Fermi liquid case M = 1 of the symmetric Anderson 
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model (2E d = —U) that the exact d-electron propagator Gd a (u) and the d- 
electron self-energy T> da (uj) = uj — Grf CT (a;) _1 obey the following local Fermi 
liquid relations in the limit w — > — iO, T — ► 0, 

Luttinger theorem: J duj G d(7 (cj) = (3.17) 

_. . , , T , 1 t ^v/n 1.9 /" 7m d\ \( W \2 /7rT\2-i 

Friedel-Langreth: -ImG dff (w) = - sm z (^— J - c |^— J + (^— J 

(3.18) 



Im£ dc7 (u;) = — T — — + c 



sm 2 (7rn d /N) V s in 2 (im d /N) > \\T K J \T K > 



(3.19) 



where c is a constant of 0(1). Combining Eqs. (3.16), (3.18) it follows 
that (for M = 1, away from particle hole symmetry) S ccr exhibits (in an 
exact theory) local Fermi liquid behavior as well, ImS co -(w — iO, T = 0) = 
a+b{u /Tk) for w — ► 0. Note that this quantity is different from the grand- 
canonical conduction electron self-energy and has a finite imaginary part 
at the Fermi level. 

The momentum dependent conduction electron Green's function in the 
presence of a single impurity is given in terms of the canonical ci-electron 
propagator as 

G Cfl(r (k,k' ; iu n ) =Gc tUT (k,iu n ) + \V\ 2 G d (iuj n )G° c ^ a (k' , iu n ) . (3.20) 

The latter expression is the starting point for treating a random system of 
many Anderson impurities [28]. 



3.2. Analytical properties and infrared behavior 
The Green's functions Gff, c have the following spectral representations 

GfteiM = <^' Af ' bAUJ l (3-21) 
with the normalization of the spectral functions A f^ jC 

/oo 
diuA fAc (uj) = 1. (3.22) 
-oo 

Taking the limit A — ► oo has important consequences on the analytical 
structure of the auxiliary particle Green's functions: 
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(1) — It follows directly from the definitions Eqs. (3.12), (3.9), using 
Eqs. (3.3), (3.4), that the traces appearing in the canonical functions Gf t b 
are taken purely over the Q = sector of Fock space 1 . Thus, the backward- 
in-time (ri < r 2 ) or hole-like contribution to the auxiliary propagators in 
Eq. (3.9) vanishes after projection, and we have 

G f<T (n - r 2 ) = -6(7i - r 2 ) km </ CT (7i)/t(T 2 )) G (3.23a) 

A^oo 

Gbfiin - r 2 ) = -G(n - r 2 ) lim (6 /i (r 1 )6t( T2 )) G . (3.23b) 

A^oo 

Consequently, their spectral functions ^4/^ have the Lehmann representa- 
tion 

A fa (u) = £ e"^ I (l,n|/t|0,m) | 2 5( w - (^ - E°J) (3.24a) 

m,ri>0 

4>» = Y, I <M|#,™> \ 2 5{u-{E l n -El)) , (3.24b) 

m,ri>0 

where £^ are the energy eigenvalues (Eq < E® is the ground state energy) 
and \Q,n) the many-body eigenstates of H in the sector Q of Fock space. 
At zero temperature, Af reduces to Af a (u) = J2 n >o I (-*-' n \fl 1 , 0) | 2 <5(lj — 
(E\ — E^)) and similar for Af,. It is seen that the Aff, have threshold 
behavior at co = Eq = E^ - Eq, with A ftb (uj) = for uj < Eq, T = 0. 
The vanishing imaginary part at frequencies u < may be shown to be a 
general property of all quantities involving slave particle operators, e.g. also 
of auxiliary particle self-energies and vertex functions. 

(2) — As will be seen in section 4.2 (Eq. (4. 2d)), physical expectation 
values not only involve the particle-like auxiliary propagators Eq. (3.23) 
but also hole-like contributions. It is, therefore, useful to define the "anti- 
fermion" and "anti-boson" propagators (in imaginary time representation) 

Gj a (n - r 2 ) = 6(r 2 - n ) lim {f^)f a (n)) G (3.25a) 

G^(n - r 2 ) = 6(r 2 - n ) lim <&tfa)6p(Ti)> G , (3.25b) 

whose spectral functions have the Lehmann representations 

A- a (u) = Y e"^ I (0,n\f a \l,m) | 2 S(u - (E° n - E^)) (3.26a) 

m,n>0 

1 This means that the auxiliary particle propagators are not calculated in the canonical 
(Q = 1) ensemble. The projection onto the Q = 1 sector of Fock space is achieved 
only when they are combined to calculate expectation values of physically observable 
operators like Gda, (S) etc. The latter can be seen explicitly, e.g., from Eq. (4. 2d), 
2nd equality. 
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^» = E I (O.nltpll.m) | 2 <5(^ - (E° n - <)) .(3.26b) 

m,n>0 

is the ground state energy in the Q = 1 sector. The expressions (3.24) 
and (3.26) immediately imply a relation between A^ b and Aj b , 

Aj tb (u) = e-^A f>b (u) . (3.27) 

(3) — The property of only forward-in-time propagation (Eqs. (3.23), 
(3.24)) means that the auxiliary particle propagators Gf^ are formally iden- 
tical to the core propagators of the well-known X-ray threshold problem 
[29]-[31]. Thus, the knowledge of the infrared behavior of the latter may 
be directly applied to the former. In particular, the spectral functions are 
found (see below) to diverge at the threshold Eq in a power law fashion 
(infrared singularity) 

A fib (u) ~ \u - E \- a f' b 9(u - Eq) (3.28) 

due to a diverging number of particle-hole excitation processes in the con- 
duction electron sea as lo — ► E Q . 

For the single channel case (M = 1), i.e. the usual Kondo or mixed 
valence problem, the exponents a/ and a b can be found analytically from 
the following chain of arguments: Anticipating that in this case the impu- 
rity spin is completely screened by the conduction electrons at temperature 
T = 0, leaving a pure-potential scattering center, the ground state |1,0) 
is a slater determinant of one-particle scattering states, characterized by 
scattering phase shifts rj a in the s-wave channel (assuming for simplicity a 
momentum independent hybridization matrix element V). To calculate the 
fermion spectral function Af a (uj) at T = from Eq. (3.24a), one needs to 
evaluate (l,n | f£ \ 0,0), which is just the overlap of two slater determi- 
nants, an eigenstate of the fully interacting Kondo system, |l,n), on the one 
hand, and the ground state of the conduction electron system in the absence 
of the impurity combined with the decoupled impurity level occupied by an 
electron with spin a, /J|0, 0), on the other hand. As shown by Anderson 
[29], the overlap of the two ground state slater determinants, (1, | fl \ 0, 0), 
tends to zero in the thermodynamic limit (orthogonality catastrophe). Anal- 
ogous relations hold for the boson spectral function A b (uj). As a result, the 
long-time relaxation into the interacting ground state is inhibited, leading 
to the infrared power law divergence of the spectral functions, Eq. (3.28). 

The X-ray threshold exponents can be expressed in terms of the scat- 
tering phase shifts at the Fermi level by the exact relation [32] 

«^ = 1-E(^) 2 - (3-29a) 
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Here the r]f a > are the scattering phase shifts of the single-particle 

wave functions in channel a' of the fully interacting ground state |1,0), 
relative to the wave functions of the free state /J|0, 0) (M|0, 0) ). Via the 
Friedel sum rule, the scattering phase shifts are, in turn, related to the 
change An ccr / of the average number of conduction electrons per scattering 
channel a 1 due to the presence of the impurity: r/f^ a > = 7rAn CCT '. Obviously, 
An ccr / is equal and opposite in sign to the difference of the average impurity 
occupation numbers of the states |1, 0) and /J|0, 0) ( 6^|0, 0) ). Thus, in the 
pseudofermion propagator Gf a we have the phase shifts, 

Vfa> = _7r (^ ~ <W) (3.30a) 

and in the slave boson propagator Gb, 

Vb = , (3.31a) 

where rid denotes the total occupation number of the impurity level in the 
interacting ground state. (The term 5 aa i in rjf a t appears because f%-\0, 0) has 
impurity occupation number 1.) For example, in the Kondo limit re^ — ► 1 
and for a spin 1/2 impurity (N = 2) this leads to resonance scattering, 
riffia' = 7r /2- As a result, one finds [33] for the threshold exponents 

a, = (3.32a) 

2 

a b = 1 - ^ (3.32b) 

These results have been found independently from Wilson's numerical renor- 
malization group approach [34, 35] and using the Bethe ansatz solution 
and boundary conformal field theory [36]. It is interesting to note that 
(i) the exponents depend on the level occupancy (in the Kondo limit 
rid ~ ¥ 1) a f = V-W an d O-b = 1 — 1/-/V, whereas in the opposite, empty 
orbital, limit rid ~ ^ 0, af — > and ab — > 1) (ii) the sum of the exponents 
a f + a b = l + 2 r -^f^ > 1. 

We stress that the above derivation of the infrared exponents a/^ holds 
true only if the impurity complex acts as a pure potential scattering center 
at T = 0. This is equivalent to the statement that the conduction electrons 
behave locally, i.e. at the impurity site, like a Fermi liquid. Conversely, 
in the multi-channel (non-FL) case, N > 2, M > N , the exponents have 
been found from a conformal field theory solution [23] of the problem in the 
Kondo limit to be a f = M/(M + N), a b = N/(M + N), which differ from 
the FL values. Thus, one may infer from the values of af^ as a function of 
rid, whether or not the system is in a local Fermi liquid state. 
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4. Mean field and non-crossing approximations 

For physical situations of interest, the s—d hybridization of the Anderson 
model (2.1) is much smaller than the conduction band width, Af(0)V <C 1, 
where A/"(0) = 1 /D is the local conduction electron density of states at the 
Fermi level. This suggests a perturbation expansion in J\f(0)V. A straight- 
forward expansion in terms of bare Green's functions is not adequate, as 
it would not allow to capture the physics of the Kondo screened state, or 
else the infrared divergencies of the auxiliary particle spectral functions 
discussed in the last section. In the framework of the slave boson represen- 
tation, two types of nonperturbative approaches have been developed. The 
first one is mean field theory for both the slave boson amplitude (b) and the 
constraint ((Q) = 1 rather than Q = 1). The second one is resummation of 
the perturbation theory to infinite order. 

4-1. Slave boson mean field theory 

Slave boson mean field theory is based on the assumption that the slave 
bosons condense at low temperatures such that (bp) / 0. Replacing the 
operator bp in H+XQ by (bp) (see Ref. [37]), where A is a Lagrange multiplier 
to be adjusted such that (Q) = 1, one arrives at a resonance level model for 
the pseudofermions. The position of the resonance, Ed + A, is found to be 
given by the Kondo temperature Tk, and is thus close to the Fermi energy. 
The resonance generates the low energy scale Tk, and leads to local Fermi 
liquid behavior. While this is qualitatively correct in the single-channel 
case, it is in blatant disagreement with the exactly known behavior in the 
multi-channel case. The mean field theory can be shown to be exact for 
M = 1 in the limit N — ► oo for a model in which the constraint is softened 
to be Q = N/2. However, for finite A?" the breaking of the local gauge 
symmetry, which would be implied by the condensation of the slave boson 
field, is forbidden by Elitzur's theorem [38]. It is known that for finite N the 
fluctuations in the phase of the complex expectation value (bp) are divergent 
and lead to the suppression of (6^) to zero (see also [39]-[41]). This is true 
in the cartesian gauge, whereas in the radial gauge the phase fluctuations 
may be shown to cancel at least in lowest order. It has not been possible to 
connect the mean field solution, an apparently reasonable description at low 
temperatures and for M = 1, to the high temperature behavior (T S> Tk), 
dominated by logarithmic temperature dependence, in a continuous way 
[37]. Therefore, it seems that the slave boson mean field solution does not 
offer a good starting point even for only a qualitatively correct description 
of quantum impurity models. 
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Fig. 2. Diagrammatic representation of the generating functional <& of the NCA. 
Also shown are the pscudoparticlc self-energies and the local electron Green's func- 
tion derived from <J>, Eqs. (19)-(21). Throughout this article, dashed, wavy and 
solid lines represent fermion, boson, and conduction electron lines, respectively. In 
the diagram for S f a the spin labels are shown explicitly to demonstrate that there 
are no coherent spin fluctuations taken into account. 

4-2. 1/N expansion vs. self-consistent formulation 

The critical judgement of mean field theory is corroborated by the results 
of a straightforward 1/iV-expansion in the single channel case, keeping the 
exact constraint, and not allowing for a finite bose field expectation value 
[42]. Within this scheme the exact behavior of the thermodynamic quanti- 
ties (known from the Bethe ansatz solution) at low temperatures as well as 
high temperatures is recovered to the considered order in 1/N. Also, the 
exact auxiliary particle exponents exf^ are reproduced in order 1/N, using 
a plausible exponentiation scheme [43]. 

In addition, dynamical quantities like the d-electron spectral function 
and transport coefficients can be calculated exactly to a desired order in 
1/N, within this approach. However, as clear-cut and economical this 
method may be, it does have serious limitations. For once, the experi- 
mentally most relevant case of N = 2 or somewhat larger is not accessible 
in 1/N expansion. Secondly, non- Fermi liquid behavior, being necessar- 
ily non-perturbative in 1/N, cannot be dealt with in a controlled way on 
the basis of a l/A^-expansion. To access these latter two regimes, a new 
approach non-perturbative in 1/N is necessary. 

We conjecture that this new approach is gauge invariant many-body the- 
ory of pseudofermions and slave bosons. As long as gauge symmetry violat- 
ing objects such as Bose field expectation values or fermion pair correlation 
functions do not appear in the theory, gauge invariance of physical quan- 
tities can be guaranteed in suitably chosen approximations by the proper 
match of pseudofermion and slave boson properties, without introducing 
an additional gauge field. This requires the use of conserving approxima- 
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tions [44], derived from a Luttinger-Ward functional <!>. consists of all 
vacuum skeleton diagrams built out of fully renormalized Green's functions 
Gbj, c and the bare vertex V. The self-energies £&,/, c are obtained by tak- 
ing the functional derivative of <3> with respect to the corresponding Green's 
function (cutting the Green's function line in each diagram in all possible 
ways) , 

£ bi/ , c = 5*/8G bJiC . (4.1) 

Irreducible vertex functions, figuring as integral kernels in two-particle Bethe- 
Salpeter equations, are generated by second order derivatives of $. 

The choice of diagrams for defines a given approximation. It should 
be dictated by the dominant physical processes and by expansion in a small 
parameter, if available. As noted before, in the present context, we may take 
the hybridization V to be a small quantity (dimensionless parameter N V). 
This suggests to start with the lowest order (in V) diagram of which is 
second order (see Fig. 2). The self-energies generated from this obey after 
analytic continuation to real frequencies (iuj — ► to — iO) and projection the 
following equations of self-consistent second order perturbation theory 



T$ CA \u> - iO) = T]T| ± [1 - f(e)]A^ a (e)G b ^ -e-iO) (4.2a) 
T$ GA \u - iO) = T W ^ f(e)A% a (e)G fa (u + e - iO) (4.2b) 

a J ™ 

G$? A \»-«>) = Jdee-^[G fa (u; + e-iO)A b[i (e) 

-A fa (e)G bil (e-cj + iO)} (4.2c) 
= Jde [G fa (u + e - iQ)A^(e) - Aj a (e)G hfl (e -u + iO)] , 

where A® = ^lmG^ a /Af(0) is the (free) conduction electron density of 
states per spin and channel, normalized to the density of states at the Fermi 
level A/"(0), and /(e) = l/(exp(/3e) + 1) denotes the Fermi distribution func- 
tion. Together with the expressions (3.10), (3.11) for the Green's functions, 
Eqs. (4.2a)-(4.2d) form a set of self-consistent equations for Sbj )C , com- 
prised of all diagrams without any crossing propagator lines and are, thus, 
known as the "non-crossing approximation", in short NCA [45, 46]. 

At zero temperature and for low frequencies Eqs. (4.2a) and (4.2b) may 
be converted into a set of linear differential equations for Gf and G b [47], 
which allow to find the infrared exponents as aj = M ^ N ; a b = M ^_ N , 
independent of n^. For the single channel case these exponents do not 
agree with the exact exponents derived in section 3. This indicates that 
the NCA is not capable of recovering the local Fermi liquid behavior for 
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M = 1. A numerical evaluation of the d-electron Green's function, which 
is given by the local self-energy S c divided by V 2 and hence is given by 
the boson-fermion bubble within NCA (Fig. 2), shows indeed a spurious 
singularity at the Fermi energy [48]. The NCA performs somewhat better 
in the multi-channel case, where the exponents a/ and at, yield the correct 
non-Fermi liquid exponents of physical quantities as known from the Bethe 
ansatz solution [19] and conformal field theory [23]. However, the specific 
heat and the residual entropy are not given correctly in NCA. Also, the 
limiting low temperature scaling laws for the thermodynamic quantities are 
attained only at temperatures substantially below Tk , in disagreement with 
the exact Bethe ansatz solution. 



4-3. Low-temperature evaluation of the self-consistency equations 



In order to enter the asymptotic power law regime of the auxiliary spec- 
tral functions the self-consistent scheme, in particular the NCA, must be 
evaluated for temperatures several orders of magnitude below Tk, the low 
temperature scale of the model. The equations are solved numerically by 
iteration. In the following we describe the two main procedures to make the 
diagrammatic auxiliary particle technique suitable for the lowest tempera- 
tures. 

The grand-canonical expectation value of the auxiliary particle number 
appearing in Eq. (3.8) is given in terms of the grand-canonical (unprojected) 
auxiliary particle spectral functions Af t b(oJ, A) by, 

(Q)g(A) = f duj [f{uj) ^Afriu, A) + 6H ]T A bfi (u, A)] (4.3) 

where f(u>), b(u) are the Fermi and Bose distribution functions, respec- 
tively. Substituting this into the expression (3.6) for the canonical partition 
function we obtain after carrying out the transformation uj — > uj + A, and 
taking the limit A — > oo 

e -m mp (T) = _p_ = lim e ^(Q) G (X) 



A^oo 



= / due-** []T A fa {u) + £ Mai?)} ■ (4-4) 



By definition Fi mp = — jj ln(Zc / Zq = q) is the impurity contribution to the 
Free energy. 

The numerical evaluation of expectation values like (Q)g(A — > oo) (Eq. 
(4.4)) or Grf MO -(o;,A — ► oo) (Eq. (4. 2d)) is non-trivial, (1) because at T = 
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the auxiliary spectral functions Aj %b {us,T) are divergent at the threshold 
frequency Eq, where the exact position of Eq is a priori not known, and (2) 
because the Boltzmann factors eT^ diverge strongly for us < 0. Therefore, 
we apply the following transformations: 

(1) Before performing the projection us ^ us + X, X ^ oo we re-define 
the frequency scale of all auxiliary particle functions A f >b according to us — > 
us+Xq, where Ao is a finite parameter. In each iteration Ao is then determined 
such that 

/ duse-P"[J2A fa (u)+J2A bfi (us)] =1 (4.5) 

(7 /I 

where Af jb (us) = lim^oo Af^(us+ Ao, A) is now an auxiliary spectral function 
with the new reference energy. It is seen by comparison with Eq. (4.4) that 
A (T) = F imp (T) = Fq =1 (T)-F q=0 (T), i.e. A is the chemical potential for 
the auxiliary particle number Q, or equivalently the impurity contribution 
to the Free energy. The difference of the Free energies becomes equal to the 
threshold energy Eq = Eq^ — Eq= at T = 0. More importantly, however, 
the above way of determining a "threshold" is less ad hoc than, for example, 
defining it by a maximum in some function appearing in the NCA equations. 
It is also seen from Eq. (4.5) that this procedure defines the frequency 
scale of the auxiliary particles such that the T = threshold divergence of 
the spectral functions is at the fixed frequency us = 0. This substantially 
increases the precision as well as the speed of numerical evaluations. Eq. 
(3.8) for the projected d-electron Green's function becomes 

G d (us) = hm eP x G d (u,T,\). (4.6) 

A^oo 

(2) The divergence of the Boltzmann factors implies that the self-consis- 
tent solutions for Af^(us) vanish exponentially ~ for negative frequen- 
cies, confirming their threshold behavior. It is convenient, not to formulate 
the self-consistent equations in terms of A like in earlier evaluations [49] , 
but to define new functions Af^(us) and ImSj^u;) such that 

A fjb {u) = f(-co) A f , b (us) (4.7) 
Im£ /;b H = f(-uj) lmt Lb (us). (4.8) 

After fixing the chemical potential Ao and performing the projection onto 
the physical subspace, the canonical partition function (Eq. (3.6)) behaves 
as limA-^oo e' 3( - A ~ Ao ^ Zc(T) = 1. In this way all exponential divergencies 
are absorbed by one single function for each particle species. The NCA 
equations in terms of these functions are free of divergencies of the statistical 
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factors and read 



lmt fa (u-i0) = Tj2j de I{ £) \_l^ ) £)) A% a (s)A m (u - e) (4.9) 
ImS^u-iO) = Tj2j ds^^j^±^Al a (e)A fa (u + e) (4.10) 
(Q)(A ,A^oo) = f dw/Hf^^H + S^H] =1 ( 4 - n ) 

ImG^w-iO) = y"de[/(e + a;)/(-e) + /(-e - w)/(e)]i /tr (e + ^(e). 

(4.12) 

The real parts of the self-energies Ej, are determined from ImSj, ImS?, 
through a Kramers-Kroenig relation, and the auxiliary functions Af a (uj) = 
^ImS /CT (w - i0)/[(w + A - i0 - £ d - ReS /CT (cj - iO)) 2 + ImS /(J (w - iO) 2 ], 
i^H = ilmS 5/i (^-fO)/[(LJ + Ao-zO-ReS5 /i (^-iO)) 2 + ImS5 /i (to-iO) 2 ], 
thus closing the above set of equations. 

The method described above allows to solve the NCA equations effec- 
tively for temperatures down to typically T = 10 -4 Tr-. It may be shown 
that the same procedure can also be applied to self-consistently compute 
vertex corrections beyond the NCA (see section 5). 

5. Conserving T-matrix approximation 

5.1. Dominant contributions at low energy 

In order to eliminate the shortcomings of the NCA mentioned above, 
the guiding principle should be to find contributions to the vertex functions 
which renormalize the auxiliary particle threshold exponents to their cor- 
rect values, since this is a necessary condition for the description of FL and 
non-FL behavior, as discussed in section 3. Furthermore, it is instructive to 
realize that in NCA any coherent spin flip and charge transfer processes are 
neglected, as can be seen explicitly from Eqs. (4.2a), (4.2b) or from Fig. 2. 
These processes are known to be responsible for the quantum coherent col- 
lective behavior of the Anderson impurity complex below Tk- The existence 
of collective excitations in general is reflected in a singular behavior of the 
corresponding two-particle vertex functions. In view of the tendency of 
Kondo systems to form a collective spin singlet state, we are here interested 
in the spin singlet channel of the pseudofermion-conduction electron vertex 
function and in the slave boson-conduction electron vertex function. It may 
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TCO t'co' fc f f f 



a u co-£2 a |i' co'-£2 c b c c c 




VCO v'co' be b b b 



Fig. 3. Diagrammatic representation of the Bethe-Salpeter equation for (1) the 
conduction electron-pseudofermion T-matrix Eq. (5.1), and (2) the conduc- 

tion electron-slave boson T-matrix T< cb \ Eq. (5.2). is obtained from T^^ 

by interchanging / <-> b and c <-> s . 

be shown by power counting arguments (compare appendix A) that there 
are no corrections to the NCA exponents in any finite order of perturbation 
theory [50]. Thus, we are led to search for singularities in the aforemen- 
tioned vertex functions arising from an infinite resummation of terms. From 
the preceding discussion it is natural to perform a partial resummation of 
those terms which, at each order in the hybridization V, contain the maxi- 
mum number of spin flip or charge fluctuation processes, respectively. This 
amounts to calculating the conduction electron-pseudofermion vertex func- 
tion in the "ladder" approximation definied in Fig. 3, where the irreducible 
vertex is given by V 2 Gt- In analogy to similar resummations for an inter- 
acting one-component Fermi system, we call the total c-f vertex function 
T-matrix T^ c f\ The Bethe-Salpeter equation for T^ c ^ reads (Fig. 3 (1)), 

T^l^iuJn, iJ n , itt n ) = + V 2 G bfl (iuj n + %J n - in n )5 aT '5 Tr7 > 

- V 2 Tj2G bfl (iuj n + iul - in n ) x (5.1) 

A similar integral equation holds for the charge fluctuation T-matrix 
(Fig. 3 (2)). 

- V 2 TJ2 G fa {iu n + iu% - m n ) x (5.2) 
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Fig. 4. Imaginary part of the conduction electron-pseudofermion T-matrix T^ c ^ 
as a function of the COM frequency O for the single-channel case M — 1, TV = 2, 
evaluated by inserting NCA solutions for the intermediate state propagators {Ed = 
-0.67D, T = 0.15D, T = 4 • lO^T^). The contribution from the pole positioned 
at a negative frequency ft — Q c f ~ — Tk (compare text) is clearly seen. 

G bjl {i^) G° cl/a (-i^-in n ) T^J v ,{iJ^iJ n ,iSl n ). 

In the above Bethe-Salpeter equations a, r, a' , r' represent spin and /i, 
^, /[/, z/ channel indices. Inserting NCA Green's functions for the interme- 
diate state propagators of Eq. (5.1) and solving it numerically, we find at 
low temperatures and in the Kondo regime (n^ ~ 0.7) a pole of in 
the singlet channel (see appendix A) as a function of the center-of-mass 
(COM) frequency f2, at a frequency which scales with the Kondo tempera- 
ture, 0, = £l c f ~ —Tk- This is shown in Fig. 4. The threshold behavior of 
the imaginary part of as a function of U with vanishing spectral weight 
at negative frequencies and temperature T = is clearly seen. In addition, 
a very sharp structure appears, whose broadening is found to vanish as the 
temperature tends to zero, indicative of a pole in T^ c ^ at the real frequency 
S7 c y , i.e. the tendency to form a collective singlet state between the conduc- 
tion electrons and the localized spin. Similarly, the corresponding T-matrix 
T( cb ) i n the conduction electron-slave boson channel, evaluated within the 
analogous approximation, develops a pole at negative values of SI in the 
empty orbital regime (n^ ~ 0.3). In the mixed valence regime (n^ ~ 0.5) 
the poles in both T^ c ^ and coexist. The appearance of poles in the 
two-particle vertex functions T^ c ^ and T^ cb \ which signals the formation 
of collective states, may be expected to influence the behavior of the system 
in a major way. 
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5.2. Self-consistent formulation: CTMA 



On the level of approximation considered so far, the description is not 
yet consistent: In the limit of zero temperature the spectral weig ht of 
and at negative frequencies f2 should be strictly zero (threshold prop- 
erty). Nonvanishing spectral weight at < like a pole contribution for 
negative in or would lead to a diverging contribution to the 

self-energy, which is unphysical. However, recall that a minimum require- 
ment on the approximation used is the conservation of gauge symmetry. 
This requirement is not met when the integral kernel of the T-matrix equa- 
tion is approximated by the NCA result. Rather, the approximation should 
be generated from a Luttinger-Ward functional. The corresponding gener- 
ating functional is shown in Fig. 5. It is defined as the infinite series of all 
vacuum skeleton diagrams which consist of a single ring of auxiliary par- 
ticle propagators, where each conduction electron line spans at most two 
hybridization vertices (Fig. 5). As shown in appendix A by means of a can- 
cellation theorem, the CTMA includes, at any given loop order, all infrared 
singular contributions to leading and subleading order in the frequency u. 
The first diagram of the infinite series of CTMA terms corresponds to NCA 
(Fig. 2). The diagram containing two boson lines is excluded since it is 
not a skeleton. Although the spirit of the present theory is different from 
a large N expansion, it should be noted that the sum of the $ diagrams 
containing up to four boson lines includes all terms of a 1/N expansion up 
to 0(1/N 2 ) [51]. By functional differentiation with respect to the conduc- 
tion electron Green's function and the pseudofermion or the slave boson 
propagator, respectively, the shown <I> functional generates the ladder ap- 
proximations T( c f\ T"( c6 ) for the total conduction electron-pseudofermion 
vertex function (Fig. 3) and for the total conduction electron-slave boson 
vertex function. The auxiliary particle self-energies are obtained in the 
conserving scheme as the functional derivatives of $ with respect to Gf or 
Gb, respectively (Eq. (4.1)). This defines a set of self-consistency equa- 
tions, which we term conserving T-matrix approximation (CTMA), where 
the self-energies are given as nonlinear and nonlocal (in time) functionals 
of the Green's functions, while the Green's functions are in turn expressed 
in terms of the self-energies. The CTMA equations are derived explicitly in 
appendix B. The solution of these equations requires that the T-matrices 
have vanishing spectral weight at negative COM frequencies Q. Indeed, the 
numerical evaluation shows that the poles of and are shifted to 
= by self-consistency, where they merge with the continuous spectral 
weight present for 17 > 0, thus renormalizing the threshold exponents of the 
auxiliary spectral functions. 

The self-consistent solutions are obtained by first solving the linear 



24 



CTMA PRINTED ON FEBRUARY 1, 2008 




+ 



spin charge 
fluctuations 

Fig. 5. Diagrammatic representation of the Luttinger-Ward functional generat- 
ing the conserving T-matrix approximation (CTMA). The terms with the con- 
duction electron lines running clockwise (labelled "spin fluctuations") generate 
the T-matrix while the terms with the conduction electron lines running 

counter-clockwise (labelled "charge fluctuations") generate the T-matrix T^ cb \ 

Bethe-Salpeter equations for the T-matrices by matrix inversion, com- 
puting the auxiliary particle self-energies from T^) and T (cb \ and then 
constructing the fermion and boson Green's functions from the respective 
self-energies. This process is iterated until convergence is reached. We 
have obtained reliable solutions down to temperatures of the order of at 
least 10~ 2 T K both for the single-channel and for the two-channel Anderson 
model. Note that Tk —> in the Kondo limit; in the mixed valence and 
empty impurity regimes, significantly lower temperatures may be reached, 
compared to the low temperature scale of the model. 

As shown in Fig. 6 (1) a), the auxiliary particle spectral functions ob- 
tained from CTMA [52] are in good agreement with the results of a nu- 
merical renormalization group (NRG) calculation [34] (zero temperature 
results), given the uncertainties in the NRG at higher frequencies. Typ- 
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Fig. 6. (1) Pseudofermion and slave boson spectral functions Af and At in the 
Kondo regime (N = 2; E d = -0.05, Y = 0.01 in units of the half bandwidth D), 
for a) the single-channel (M = 1) and b) the multi-channel (M = 2) case. In 
a) the symbols represent the results of NRG for the same parameter set, T = 0. 
The slopes of the dashed lines indicate the exact threshold exponents as derived in 
section 3 for M = 1 and as given by conformal field theory for M = 2. The insets 
show magnified power law regions. (2) CTMA results (symbols with error bars) 
for the threshold exponents a/ and ab of Af and Af,, N = 2, M = 1. Solid lines: 
exact values (section 3), dashed lines: NCA results (section 4.2). 



ical behavior in the Kondo regime is obtained: a broadened peak in Af, 
at lo ~ | -Ed |, representing the hybridizing d-level and a structure in Af at 
uj ~ Tk- Both functions display power law behavior at frequencies below 
Tk, which at finite T is cut off at the scale to ~ T. The exponents extracted 
from the frequency range T < co < Tk of our finite T results compare 
well with the exact result also shown (see insets of Fig. 6 (la)). A similar 
analysis has been performed for a number of parameter sets spanning the 
complete range of d-level occupation numbers rid- The extracted power law 
exponents are shown in Fig. 6 (2), together with error bars estimated from 
the finite frequency ranges over which the fit was made. The comparatively 
large error bars in the mixed valence regime arise because here spin flip 
and charge fluctuation processes, described by the poles in and T^ cb \ 
respectively, are of equal importance, impeding the convergence of the nu- 
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Fig. 7. Static susceptibility of the two-channel Anderson impurity model: CTMA 
and NCA results (E d = -0.8D, T = 0.1D, Landc factor g = 2). 



merical procedure. In this light, the agreement with the exact results (solid 
curves) is very good, the exact value lying within the error bars or very 
close in each case. 

In the multi-channel case (N > 2, M > N) NCA has been shown 
[50] to reproduce asymptotically the correct threshold exponents, aj = 
M/(M + N), a b = N/(M + N), in the Kondo limit. Calculating the T- 
matrices using NCA Green's functions (as discussed in the single-channel 
case) we find again a pole in the singlet channel of However, in this 

case the CTMA does not renormalize the NCA exponents in the Kondo 
limit of the two-channel model, i.e. the threshold exponents obtained from 
the CTMA solutions are very close to the exact ones, aj = 1/2, a b = 1/2, 
as shown in Fig. 6 (1) b). 

The agreement of the CTMA exponents with their exact values in the 
Kondo, mixed valence and empty impurity regimes of the single-channel 
model and in the Kondo regime of the two-channel model may be taken 
as evidence that the T-matrix approximation correctly describes both the 
FL and the non-FL regimes of the SU(N)xSU(M) Anderson model (N=2, 
M=l,2). Therefore, we expect the CTMA to correctly describe physi- 
cally observable quantities of the SU(N) x SU(M) Anderson impurity model 
as well. As a check, we have calculated the static spin susceptibility x 
of the two-channel Anderson model in the Kondo regime by solving the 
CTMA equations (see appendix B) in a finite magnetic field H coupled 
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to the impurity spin and taking the derivative of the magnetization M = 
5<7MB( n /T — nfi) with respect to H. The resulting x(T) = (dM/dH)x is 
shown in Fig. 7. It is seen that CTMA correctly reproduces the exact [19] 
logarithmic temperature dependence below the Kondo scale Tk- In con- 
trast, the NCA solution recovers the logarithmic behavior only far below 
Tk- Other physical quantities will be calculated for the Anderson model in 
forthcoming work. 

6. Conclusion 

We have presented a novel technique to describe correlated quantum im- 
purity systems with strong onsite repulsion, which is based on a conserving 
formulation of the auxiliary boson method. The conserving scheme allows 
to preserve the conservation of the local charge Q without taking into ac- 
count time dependent fluctuations of the gauge field A. Taking, as a result, 
A as a time independent quantity represents a great simplification of this 
approach. As a standard diagram technique this method has the potential 
to be applicable to problems of correlated systems on a lattice as well, while 
keeping the full dynamics of the pseudofermion and slave boson fields. By 
including the leading infrared singular contributions (spin flip and charge 
fluctuation processes), this technique allows to correctly describe physical 
quantities, like the magnetic susceptibility, both in the Fermi and in the 
non-Fermi liquid regime, over the complete temperature range, including 
the crossover to the correlated many-body state at the lowest temperatures. 

We wish to thank S. Bocker, T.A. Costi, S. Kirchner, A. Rosch, A. Ruck- 
enstein and Th. Schauerte for stimulating discussions. S. Bocker has per- 
formed part of the numerical solutions. This work is supported by DFG 
through SFB 195 and by the Hochleistungsrechenzentrum Jiilich through a 
grant of computer time on a Cray T3E parallel computer. 



Appendix 

A. Infrared cancellation of non CTMA diagrams 

The CTMA is not only justified on physical grounds by the inclusion of 
the maximum number of spin flip and charge fluctuation processes at any 
given order of perturbation theory, but also by an infrared cancellation of 
all diagrams not included in the CTMA. In the following we will prove this 
cancellation theorem. 

(1) Power counting. — Each auxiliary particle loop carries a factor of 
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Fig. 8. Typical pseudofermion self-energy skeleton diagram of loop order L = 6, 
containing L sp = 1 spin (or fermion) loop and L c h = 1 channel loop. 



the fugacity exp(— f3X), which vanishes upon projection onto the Q = 1 
subspace, A — > oo. Therefore, an arbitrary / or b self-energy diagram 
consists of one single line of alternating fermion and boson propagators, 
with the hybridization vertices connected by conduction electron lines in 
any possible way, as shown in Fig. 8 (see also [45]). Such a fermion self- 
energy skeleton diagram of loop order L is calculated as 

yj(. L) ( w ) = (_1)3L-1+L sp N L sp M L ch pi 

x f^l...±L /( £l ) . . . f {£L ) A^e,) . . . A° c (s L s L ) (A.l) 
J vr 7T 

X G b (uJ + UJi)Gf(uJ + UJi + Uj'i) ... 

k k-l k 

x ... G b (u + ^2ui + ^2(4j Gf(u + ^2(ui + (4)) ... G b {oj + uj L ) , 

i=l i=l i=l 

where Gf tb are the renormalized, i.e. power law divergent auxiliary particle 
propagators, and L sp and L c /j denote the number of spin (or fermion, c-f) 
loops and the number of channel (or c-b) loops contained in the diagram, 
respectively. Spin and channel indices are not shown for simplicity. Each of 
the auxiliary particle frequencies cjj, u)[ coincides with one of the integration 
variables £j, j = 1, . . . ,L, in such a way that energy is conserved at each 
hybridization vertex. This implies that the sign of the frequency carried 
by a c-electron line is Sj = +, if the c-electron line runs from right to 
left, and Sj = — , if it runs from left to right in Fig. 8. An analogous 
expression holds for the slave boson self-energy diagrams. By substituting 
Xj = Sj/u, j = 1, . . . , L and factoring out uj~ a f'( bS > from each fermion (boson) 
propagator, the infrared behavior of the term Eq. (A.l) is deduced as 

Im£$M = Cu a f' b+L ^- a f- a >>\ (A.2) 

where C is a finite constant. Clearly, when the NCA solutions are inserted 
for the propagators G/,6, i.e. a/ + a b = 1, their power law behavior is just 
reproduced by any term of the form Eq. (A.l). However, this is no longer 
the case for the exact propagators in the Fermi liquid regime (M < N), 



ctma printed on February 1, 2008 



29 




Fig. 9. Set of contributions to skeleton diagrams not contained in CTMA which 
cancel in the infrared limit to leading and subleading order in the external fre- 
quency, bJ — > 0. 

where in general a/ + at, > 1. Thus, the infinite resummation of terms to 
arbitrary loop order is unavoidable in this case. 

(2) Infrared cancellation. — As discussed in section 5.2, the CTMA is 
equivalent to the self-consistent summation of all skeleton free energy di- 
agrams, where a conduction electron line spans at most two hybridization 
vertices (Fig. 5). Thus, any skeleton self-energy diagram not included in 
CTMA contains at least one conduction electron "arch" which spans four 
(or more) vertices, with four conduction lines reaching from inside to out- 
side of the arch as shown in Fig. 9(a). For each such diagram there exists 
another skeleton, which differs from Fig. 9(a) only in that the end points of 
two conduction lines inside the arch are interchanged (Fig. 9(b)). The corre- 
sponding permutation of fermionic operators implies a relative sign between 
the terms Fig. 9(a) and (b). Without loss of generality we now assume to > 
for the external frequency of the self-energy. The leading infrared singular 
behavior of the term (A.l) arises from those parts of the integrations, where 
the arguments of the Gf,Gb are such that the divergences of all propaga- 
tors lie within the integration range. This implies at least —uj < Ej < 0, 
j = 1,...,L. Therefore, the terms corresponding to Fig. 9(a), (b) differ 
only in the frequency arguments of the Green's functions inside the arch, 
and at temperature T = the leading infrared behavior of their sum reads, 

£ ( f L,o) (w) + £ ( /" 6) H W ~ (-l) 3L - 1+Lsp N Ls p M Lch T L x (A.3) 
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r dei d£L [Gf{Q + J + wi)G b (Q + J + wi + w 2 ) 

—Gf(u> + u/ + u^G^u; + u/ + u; 3 + w 2 )] 

Here u) denotes the sum of all frequencies w, £j entering the diagrammatic 
part, Fig. 9, from the left, and F(u,{£j}) consists of all terms which are 
not altered by interchanging the c-electron lines. In the infrared limit, 
ll>i — — > 0, the term in square brackets may be written as 

-^[G f {uj + u')G b {u + uj' + u 2 )]{uJi - co 3 ) (A.4) 

and upon performing the integrations over ui, UJ3 the difference {u\ — L03) 
leads to an additional factor of w. A similar cancellation of the leading 
infrared singularity occurs between the terms shown in Fig. 9(c), (d). In 
an analogous way it may be shown that combining the terms Fig. 9 (a)-(d) 
leads to a factor of lo 2 compared to the power counting result for one single 
term. Thus, the infrared singularity of all non-CTMA terms of loop order 
L is weaker than the Lth order CTMA terms by at least O(oj 2 ), 

45 a V) + • • • + 4 L,< V) ^° iv a f' b+L{1 - a f- a » )+2 . (A.5) 

It should be emphasized that in the above derivation, L appears only as a 
parameter and, thus, the cancellation theorem holds for arbitrarily high loop 
order L. This proves that the CTMA captures the leading and subleading 
infrared singularities [u — ► 0) at any given order L. 

B. CTMA equations 

In this appendix we give explicitly the self-consistent equations which 
determine the auxiliary particle self-energies within CTMA. For that pur- 
pose, it is useful to define conduction electron-fermion and conduction 
electron-boson vertex functions T^W, without (+) or with (— ) 

an alternating sign between terms with even and odd number of rungs (com- 
pare Fig. 3). In the Matsubara representation, these vertex functions, to 
be labelled "even" (+) and "odd" (— ) below, are given by the following 
Bethe-Salpeter equations: 

T (cf)(±) £(iu; n ,iu;' n ,in n ) = U {cf) £(iu n ,iu;' n ,in n ) 

±V2 \ E G bfi(i"n + iUn ~ *«n) X (B.la) 
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and 



U {cf) £{iu n , iu' n , ifl n ) = -V 4 - G b-^n + iul - in n ) x (B.lb) 

to" 

G /ff «) G% T (i£l n - O + iu£ - *n n ) 

±y2 \ E G/<x(^„ + ^ - in n ) x (B.2a) 

Gbfi{iu'n) G° cua {iuj" - ifl n ) T {cb)i± ^(iuj",iuj' n ,ittn) 
U {cb ll{iu n ,iJ n ,iSl n ) = -V^^Gf^n + iUn-Mn) x (B.2b) 

g^M G° cua (i^ - in n ) G fa {iJ n + iJ^ - in n ). 

Note that, in addition to the alternating sign, these vertex functions differ 
from the T-matrices defined in Eqs. (5.1), (5.2) in that they contain only- 
terms with two or more rungs, since the inhomogeneous parts U^ c ^ and 
lj(cb) re p resen t terms with two bosonic or fermionic rungs, respectively. The 
terms with a single rung correspond to the NCA diagrams and are evaluated 
separately (see below). 

The spin degrees of freedom of T^ c ^ W are uniquely determined by 
the spin indices a, r of the ingoing conduction electron and pseudofermion 
lines (Fig. 3). It is instructive to note that in the spin S = 1/2 case 
(N = 2) the singlet and triplet vertex functions (which correspond to the 
two-particle Green's functions in the singlet channel, (j) s ~ ^ a {T{{c a f- a — 
c- (7 f (7 )(c] r f_ a — cLo- /<£)}), and in the triplet channel with magnetic quan- 
turn number m = 0, ±1, <^ =0 ~ Y,AT{M- a + c_ CT / CT )(4 f_ c + cl CT /t)}), 
<f>m=±i ~ (T{c ± i / ± icj, i f^_i }), respectively) may be identified in the fol- 

2 2 2 1 

lowing way, 

T (cf)s = Y^Tl cf J}- ] (B.3a) 

a 

TSJ = E T t^ (+) ( B - 3b ) 

a 

rp{cf)t _ rp(cf)( + ) , u „\ 

2 > 2 

Analogous relations hold for the conduction electron-boson vertex function 
in terms of the channel degrees of freedom [/,, v. The total CTMA pseu- 
doparticle self-energies, as derived by functional differentiation from the 
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Fig. 10. Diagrammatic representation of the CTMA expressions for pscudoparticlc 
self-energies S/ and Ef,. The first term drawn on the righthand side of £/ and 

E b , respectively, is the NCA diagram. The diagrammatic parts t[°P , are 
explained in the text. 

generating functional Fig. 5, are shown in Fig. 10 and consist of three 
terms each, 

E /<r (iw n ) = V% CA \iu n ) + XfPiiUn) + S^(^n) (B.4) 

The first term of £/ and represents the NCA self-energies, Eqs. (4.2a), 
(4.2b). The second and third terms arise from the spin and the charge 
fluctuations, respectively, and are given for pseudofermions by 

Z { fP(iu n ) = M^ G° c (in n - iuj n )T[ cf \iuj n , iu n , ifl n ) (B.6a) 

P Q 

4?(^") = " M E G c(^n " to/JGiKj X (B.6b) 

and for slave bosons by 

4 C /W = F 2 -! 2 - iu n )G f (iJ n ) x (B.7a) 

T 2 (c/) «, < + < - *u, n )G°« - iu^G/KO 

4f On) = ^ 4 E G c(^n " ifin)^ (««;„, iw n , ifi n ) , (B.7b) 
P o 
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where the vertex functions appearing in these expressions are defined as 

jf /) = H+l T (cf) (+) + ?Lzl T {cf) (-) _ N uicf) (B 8a) 

T ( cf) = N + l t(c/) (+) _ N - 1 T(c/) ( _ } _ ^ (c/) ^ g ^ 



r ( cb) = M-l T(cb) (+) + M + l r(c6) ( _ } _ M u(cb) (B 9a) 



(cb) = M-l (c6) (+) _ M + l (c6) ( _ } _ (cb) 

2 o o 



(B.9b) 



These combinations of the even and odd vertex functions ensure the proper 
spin and channel summations in the self-energies. For the sake of clarity, the 
spin and channel indices as well as the frequency variables are not shown ex- 
plicitly In Eqs. (B.8a), (B.9a) the terms with two rungs, NU (cf \ MU (cb \ 
have been subtracted, since they would generate non-skeleton self-energy 
diagrams. Likewise, in Eqs. (B.8b), (B.9b) the two-rung terms have been 
subtracted in order to avoid a double counting of terms in the self-energies. 

We now turn to the analytic continuation to real frequencies of the 
expressions derived above. Transforming the Matsubara summations into 
contour integrals shows that integrations along branch cuts of auxiliary par- 
ticle Green's functions carry an additional factor exp(— (3\) as compared to 
integrations along branch cuts of physical Green's functions, which vanishes 
upon projection onto the the physical Fock space, A — > oo. Thus, as a gen- 
eral rule, only integrations along branch cuts of the c-electron propagators 
contribute to the auxiliary particle self-energies. Therefore, by performing 
the analytic continuation, itu n — > to — iO = uj in all frequency variables, we 
obtain the advanced pseudofermion self-energy, 

Sjf/V) = M J ^ f{e-uj) A° c (e - w) vrAA(0)T 1 {c/) \u,w,e) (B.lOa) 

E<£>(w) = -MT J ^ j^f(e-co) f(e' - u) x (B.lOb) 

A° c (lo - e)G b {e) ^AA(0)T 2 (cb) (e, e', e + e' - u)A c {u - e')G b (e') 
and the advanced slave boson self-energy, 

4f(u) =-NTj^J^f(s-u) f(e' - u) x (B.lla) 

A° c (e - u)G f {e) vrAA(0)r 2 (c/) ( e ,e', e + e> - u,)A° c (e> - u)G f {e>) 
S[f ( W ) = -N J^f(e-u)A c (u;-e)TrAf(0)Tl cb \u;,u;,e), (B.llb) 
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where the vertex functions are given by Eqs. (B.8), (B.9) with 

±(-r) f-f(e-n) x (B.12a) 

J 7T 

Gb . {uj + e _ f2) G/tr ( £ )^ T (n - £ )T (c/) (± ^(e, </, n) 
vrAA(O) £/ (c/)(± ^ T (u;,u/,n) = +r 2 J^- /(e-fi) x (B.12b) 
GV(w + e - f2)G /<r (e)^ T (n - e)G bfl (J + e - fi) 

and 

±(+r) [— f(e-n) x (B.13a) 
J vr 

G/ff(w + e _ n) G6p ( e )^( e - n) r (c6) a/, n) 

ttAA(O) U (cb) {± ^(co, u', fi) = -r 2 /" ^ /(e - 0) x (B.13b) 

G /ff ( W + e - f2)G 6p ( e )A^(e - 0) G /(T (u/ + s - fi) 

In the above expressions, like in the NCA equations (4.2), we have used the 
dimensionless conduction electron spectral density, A®(uj) = ^ ImG^ .(a; — 
iO)/ftf(Q), and we have suppressed obvious spin and channel indeces. All 
frequency variables are to be understood as the limit oj = u) — iO. 
The equations (B.10), (B.ll), supplemented by the vertex functions Eqs. 
(B.8), (B.9), (B.12), (B.13) form, together with the NCA contributions 
Eqs. (4.2a), (4.2b) and the definitions of the auxiliary particle Green's 
functions, Eqs. (3.10), (3.11), the closed set of self-consistent CTMA equa- 
tions [52]. It is seen that in these equations only those branches of the 
T-matrix vertex functions appear which are advanced with respect to all 
three frequency variables, although in general the T-matrix consists of 2 3 
independent analytical branches. This simplification is a consequence of 
the exact projection onto the physical sector of Fock space. Inspection of 
the analytically continued CTMA equations also shows that the slave bo- 
son self-energy is obtained from the pseudofermion self-energy, including 
the proper signs, by simply replacing Gf^Gb and inverting the frequency 
argument of A® in all expressions. 

Eqs. (B.10)-(B.13) may be rewritten in terms of the spectral functions 
without threshold, ^4/,&, in a straight-forward way as explained in section 
4.3, thus avoiding divergent statistical factors in the d-electron Green's 
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function. The CTMA equations are solved numerically by iteration. In 
each iteration step the NCA and the T-matrix contributions are computed 
separetely. We solve the linear T-matrix Bethe-Salpeter-Equations (B.12), 
(B.12) by discretization of the frequency integrals and subsequent solution 
of the resulting multidimensional linear equations. Because of the formula- 
tion given above, where first the vertex parts T^ c ^ W and T^ c ^ W , com- 
prised of all T-matrix diagrams with two or more rungs, are calculated and 
in a second step the non-skeleton two-rung contributions are subtracted 
(Eqs. (B.8), (B.9)), each term of the T-matrix equations involves at most 
one frequency integration. The resulting numerical effort is manageable: 
One complete iteration within the self-consistent scheme has been done on 
a parallel computer within approximately 5-10 s CPU time. 

REFERENCES 

[1] L.D. Landau, Sov.Phys. JETP 3, 920 (1956); 5, 101 (1957); 8, 70 (1959). 
[2] A. J. Millis, Phys. Rev. B 48, 7183 (1993). 

[3] T. Moriya and T. Takimoto, J. Phys. Soc. Japan 64, 960 (1995) 
[4] A. Rosch et al., Phys. Rev. Lett. 79, 159 (1997). 

[5] For non-FL behavior near a quantum spin-glass transition, see S. Sachdev, 
N. Read and R. Oppcrmann, Phys. Rev. B 52, 10286 (1995). 

[6] H. v. Lohneysen et al., Phys. Rev. Lett. 72, 3262 (1994); B. Bogenberger and 
H. v. Lohneysen, Phys. Rev. Lett. 74, 1016 (1995). 

[7] H. v. Lohneysen, J. Phys. Cond. Mat. 8 9689 (1996). 

[8] F. Steglich, J. Phys. Cond. Mat. 8 9909 (1996). 

[9] F. M. Grosche, S. R. Julian, N. D. Mathur and G. G. Lonzarich, Physica B 
223+224, 50 (1996). 

[10] S. R. Julian et al., J. Phys. Cond. Mat. 8 9675 (1996). 

[11] P. Nozicres and A. Blandin, J. Phys. (Paris) 41, 193 (1980). 

[12] A comprehensive overview appears in D. L. Cox and F. Zawadowski, Adv. 
Phys., in press (1998); (cond-mat/9704103). 

[13] D.L. Cox, Phys. Rev. Lett. 59, 1240 (1987); Physica C 153, 1642 (1988). 

[14] M. B. Maple et al., J. Low Temp. Phys. 95, 225 (1994); 99, 223 (1995); 
M. B. Maple et al., J. Phys. Cond. Mat. 8 9773 (1996). 

[15] R. Chain and M. B. Maple, J. Phys. Cond. Mat. 8 9939 (1996). 

[16] J. R. Schricffer and P. A. Wolff, Phys. Rev. B 149, 491 (1966). 

[17] K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 



36 



CTMA PRINTED ON FEBRUARY 1, 2008 



i] N. Andrei, K. Furuya, J.H. Lowenstein, Rev.Mod.Phys. 55, 331 (1983). 

3] P.B. Wiegmann, A.M. Tsvclik, Pis'ma Zh. Eksp. Teor. Fiz. 38, 489 (1983) 
[JETP Lett. 38, 591 (1983)]. 

)] P. Nozieres, J. Low. Temp. Phys. 17, 31 (1974). 

L] N. Andrei, C. Destri, Phys. Rev. Lett. 52, 364, (1984). 

2] A. M. Tsvelik, J. Phys. C 18 159 (1985). 

3] I. Affleck and A.W.W. Ludwig, Nucl. Phys. 352, 849, (1991); B360, 641, 
(1991); Phys. Rev. B 48, 7297 (1993). 

1] S. E. Barnes, J. Phys. F6, 1375 (1976); F7, 2637 (1977). 

5] A. A. Abrikosov, Physics 2, 21 (1965). 

3] P. Coleman, Phys. Rev. B29, 3035 (1984). 

7] K. Yamada, Prog. Theor. Phys. 53, 970 (1975); ibid. 54, 316 (1975); ibid. 55, 
1345 (1976); K. Yosida and K. Yamada, Prog. Theor. Phys. 53, 1286 (1975). 

3] D. L. Cox, N. E. Bickers and J. W. Wilkins, Phys. Rev. B 36, 2036 (1987). 

3] P. W. Anderson, Phys. Rev. Lett. 18, 1049 (1967). 

)] P. Nozieres and C. T. Dc Dominicis, Phys. Rev. 178, 1073; 1084; 1097 (1969). 

L] G. D. Mahan, Many-Particle Physics, 2nd ed., pp. 732 (Plenum Press, New 
York, 1990) gives an overview. 

2] K. D. Schottc and U. Schottc, Phys. Rev. 185, 509 (1969). 

3] B. Menge and E. Muller-Hartmann, Z. Phys. B73, 225 (1988). 

1] T.A. Costi, P. Schmittcckcrt, J. Kroha and P. Woffle, Phys. Rev. Lett. 73, 
1275 (1994). 

5] T.A. Costi, P. Schmitteckert, J. Kroha and P. Woffle, Physica (Amsterdam) 
235-240C, 2287 (1994). 

3] S. Fujimoto, N. Kawakami and S.K. Yang, J.Phys.Korea 29, S136 (1996). 

7] D. M. Newns and N. Read, J. Phys. C16, 3273 (1983); Adv. Phys. 36, 799 
(1988). 

3] S. Elitzur, Phys. Rev. D 12,3978 (1975). 

3] A. Jcvicki, Phys. Lett. B 71, 327 (1977). 

)] F. David, Comm. Math. Phys. 81, 149 (1981). 

L] I.D. Lawrie, J. Phys. A 18, 1141 (1985). 

2] B. Jin and Y. Kuroda, J. Phys. Soc. Japan 57, 1687 (1988). 

i] T. Matsuura et al., J. Phys. Soc. Japan 66, 1245 (1997). 

1] G. Baym and L.P. Kadanoff, Phys. Rev. 124, 287 (1961); G. Baym, Phys. Rev. 
127 1391 (1962). 

[45] H. Kcitcr and J. C. Kimball, J. Appl. Phys. 42, 1460 (1971); N. Grewe and 
H. Keiter, Phys. Rev B 24, 4420 (1981). 



ctma printed on February 1, 2008 



37 



[46] Y. Kuramoto, Z. Phys. B 53, 37 (1983); H. Kojima, Y. Kuramoto and M 
Tachiki, ibid. 54, 293 (1984); Y. Kuramoto and H. Kojima, ibid. 57, 95, 
(1984); Y. Kuramoto, ibid. 65, 29, (1986). 

[47] E. Miiller-Hartmann, Z. Phys. B57, 281 (1984). 

[48] T.A. Costi, J. Kroha, and P. Wolfle, PRB 53, 1850 (1996). 

[49] N. E. Bickers, Rev. Mod. Phys. 59, 845 (1987); N. E. Bickers, D. L. Cox & J. 
W. Wilkins, Phys. Rev. B36, 2036 (1987). 

[50] D. L. Cox and A. E. Ruckcnstcin, Phys. Rev. Lett. 71, 1613 (1993). 

[51] F. Anders and N. Crewe, Europhys. Lett. 26, 551 (1994); F. Anders, 
J. Phys. Cond. Mat. 7, 2801 (1995). 

[52] J. Kroha, P. Wolfle and T. A. Costi, Phys. Rev. Lett. 79, 261 (1997). 



